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We have computed the light and strange quark masses and the renormalization constants of the quark bilinear 
operators, by studying the large-p^ behaviour of the lattice quark propagator and 3-point functions. The calcula- 
tion is non-perturbatively improved, at 0{a), in the chiral limit. The method used to compute the quark masses 
has never been applied so far, and it does not require an explicit determination of the quark mass renormalization 
constant. 
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1. INTRODUCTION 

In this talk we present the results of two dif- 
ferent lattice calculations. The first one is a 
non-perturbative determination of the light and 
strange quark masses, based on the study of 
the large-p^ behaviour of the renormalized quark 
propagator. This method has never been applied 
so far, and it has the advantage of not requiring 
an explicit evaluation of the quark mass renor- 
malization constant. The second calculation is 
a non-perturbative determination of the renor- 
malization constants of the bilinear quark oper- 
ators, performed by using the non-perturbative 
renormalization method (NPM) . In both cal- 
culations we have used a sample of 100 gauge 
fields configurations generated, in the quenched 
approximation, on a 24"^ x 64 lattice, at f3 = 6.2. 
The quark propagators have been computed by 
using the non-perturbatively C'(a)-improved Al- 
pha action Q , at four different values of the light 
quark masses. 

The NPM, used to compute the lattice renor- 
malization constants, is standard by-now, and it 
will not be discussed here in further details. The 
procedure followed to determine the light quark 
masses is based on an axial- vector Ward identity, 
at zero momentum transfer, which relates the 
quark propagator, S{p), to the amputated Green 
function of the pseudoscalar density, A5 (p) , com- 



puted between external (off-shell) quark states of 
momentum p: 

2mq{fi)A5{p; = j5S{p;ny^ + 5(p;Ai)"^75 (1) 

All quantities in eq. (|l|) are renormalized in the 
same scheme, and at the same renormalization 
scale /i. In the RI-MOM scheme, the Green func- 
tion A5{p;fi) satisfies, in a fixed gauge, the fol- 
lowing condition: (1/12) Tr[75A5(p; yLt)]p2^^2 = 1. 
By tracing both sides of eq. (p with 75 and choos- 



ing = fJ-^, one obtains: 



(.) = ^Tr 



(2) 
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Thus, the renormalized quark mass can be com- 
puted by studying the large-p^ behaviour of the 
renormalized quark propagator. 

The results presented in this talk aim to be 
non-perturbatively C'(a)-improved. However, we 
have been able to improve the quark propaga- 
tor and the 3-point functions only in the chi- 
ral limit. The difficulties arise from the fact 
that these quantities are both off-shell and gauge 
non- invariant. In this case, in order to achieve 
the improvement, three distinct classes of higher- 
dimensional operators have to be included in the 
mixing with the lower-dimensional operators . 
These are: i) gauge-invariant operators, which do 
not vanish by the equation of motion; ii) gauge 
non-invariant operators, which are, however, al- 
lowed by BRST symmetry; iii) operators vanish- 
ing by the equation of motion, which can be either 
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gauge-invariant or non-invariant, and which give 
only rise to contact terms in the correlation func- 
tions. It may be shown that the off-shell 0{a)- 
improvement of the lattice action does not require 
effectively the inclusion of additional operators, 
besides the standard Clover term ||]. Thus, in the 
following, we will only discuss the improvement of 
the quark and the bilinear quark operators. 

2. QUARK MASSES 

The improvement of the lattice quark field, 
Ql, involves the mixing with both a gauge non- 
invariant operator and an operator vanishing by 
the equation of motion |^ : 

g = [l + ac'^{p + mQ) + ac^Gi^] Ql (3) 

with Zq = Zq {1 + hqara). Therefore, in order to 
improve the quark propagator, three unknown co- 
efficients, namely Zq, c'q and cngi, rnust be deter- 
mined. In principle, this could be done by study- 
ing the large-p^ behaviour of the lattice quark 
propagator. However, in practice, one finds that, 
up to very small logarithmic and power correc- 
tions, only the following combinations of coeffi- 
cients can be determined at large p^: 

(1/12) Tr [~ii>SL{v)] ^Zq[l^ 2a c^gi m] (4) 

(1/12) Tr [Sl{p)] ^ -2a c; + 2acNGi Zq (5) 

where Zq, in eqs. (^) and (||), is defined in 
the RI-MOM scheme. Therefore, in order to 
compute the renormalized, improved propagator, 
S{p), we have followed an approximate procedure 
(see ref. [Q for more details). We have subtracted, 
from the scalar part of the lattice propagator, the 
constant C'(a)-term to which it reduces, asymp- 
totically, at large p^, cfr. eq. (^. The propa- 
gator has been then renormalized by using, as 
an approximate value of the renormalization con- 
stant, the right-hand side of eq. (^). This pro- 
cedure leaves unsubtracted ©(amcNci) terms in 
the renormalized propagator, thus affecting the 
determination of the quark masses by 0{g^am) 
systematic errors. These errors, however, are ex- 
pected to be negligible in the case of light quarks. 

Once the improved renormalized quark propa- 
gator has been computed, the procedure to de- 
termine the light quark masses, from eq. (|^), is 



straightforward. The renormalized quark mass 
is extracted, as a function of the bare quark 
mass, by fitting the trace of S{p)~^ at large p^. 
This is shown in Figure |l|. Then, we extrapolate 
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Figure 1. Trace of the renormalized inverse quark 
propagator as a function of the scale. The differ- 
ent curves correspond to different values of the 
bare quark mass, the lowest one being the extrap- 
olation to the chiral limit. The solid lines repre- 
sent the values of the renormalized quark masses, 
obtained by fitting the propagator, in the large-p^ 
region, as a function of the scale, according to the 
predictions of the renormalization group equation. 

to the physical values of the light and strange 
quark masses, by using the experimental values 
of both the pseudoscalar and the vector meson 
masses, closely following the procedure described 
in ref. jsj . For the light and strange quark masses, 
in the RI-MOM scheme, at 2 GeV, we find: 

,nf I = 5.7(1.0) MeV , = 133(18) MeV (6) 

These values have been obtained in a complete 
non-perturbative way. They can be translated 
into the masses in the MS scheme, by using con- 
tinuum perturbation theory at the N^LO [^. At 
the scale fj, — 2 GeV, we find: 

TO/ = 4.6(8) MeV , m, = 108(15) MeV (7) 

3. RENORMALIZATION CONSTANTS 

The non-perturbative calculation of the renor- 
malization constants, by using the NPM Hi, is be- 
coming, at present, a standard practice in lattice 
calculations. Therefore, in this section, we limit 
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Table 1 



Values of the lattice renormalization constants obtained by using the NPM and one-loop boosted perturba- 
tion theory (BPT). The scheme- dependent renormalization constants are given in the RI-MOM scheme. 





Zv 


Za 


Zp/Zs 


Zp(2GeV) 


Z5(2GeV) 


ZT(2GeV) 


NPM 


0.757(4) 


0.794(6) 


0.78(2) 


0.44(2) 


0.56(2) 


0.87(2) 


BPT 


0.846 


0.862 


0.954 


0.587 


0.618 


0.942 



ourselves to present our results. The only point, 
which is worth to be mentioned, concerns the 
off-shell C'(a)-improvement of the bilinear quark 
operators. As discussed in ref. the operator 
Or = qPg also mixes, at Oia), with an operator 
which vanishes by the equation of motion: 



0'Y = q T{p + mo) + {-p + mo)r 



(8) 



However, this mixing, as well as the C'(a)-mixing 
of the quark field operator, does not contribute 
to the determination of the renormalization con- 
stant in the chiral limit [^,0 . 

In Table |l| we present our results for the renor- 
malization constant Zy, Za, Zg, Zp and Z^, and 
compare them with the predictions of one-loop 
(boosted) perturbation theory. Two remarks are 
in order at this point. The first one is that pre- 
vious results for the vector and axial-vector cur- 
rent renormalization constants have been also ob- 
tained, for the same action, by studying the lat- 
tice chiral Ward identities, in ref. They found 
Zv = 0.793 and Za = 0.809, in good agreement 
with our determination. The second remark con- 
cerns the sizeable discrepancy between the non- 
perturbative determination of Zp (but see also 
Zs) and the prediction of one- loop perturbation 
theory. It has been suggested that this discrep- 
ancy may be (mostly) due to the presence of a 
large, non-perturbative, power correction in the 
correlation function Ap, which has not been sub- 
tracted in order to compute Zp However, 
we find this not to be the case. To illustrate the 
point, we show in Figure the behaviour of the 
(projected) correlation functions Fp and F5, and 
the ratio Tp/Ts, as a function of the scale. No- 
tice, in particular, that in the region a^p^ <^ 1, 
which has been considered to compute the renor- 
malization constants, the ratio Tp/Ts exhibits a 
nice plateau, as indeed expected only under the 
assumption of negligible power corrections. 



♦ r,=i/z, 

■ r/r,=z,/Zp 

- Perturbat. Term 

- P.T. + Power Corr. 




Figure 2. Green functions of the pseudoscalar 
and scalar densities, and the ratio Tp/Ts, as a 
function of the scale. The dashed lines show the 
scale dependence predicted by the renormalization 
group equations at the NLO. The solid lines are 
the predictions obtained by including the effect of 
the leading power correction l/p^). 
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